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The effective interaction energy of a colloidal sphere in 
a suspension containing small amounts of non-ionic poly- 
mers and a flat glass surface has been measured and calcu- 
lated using total internal reflection microscopy (TIRM) and 
a novel approach within density functional theory (DFT), re- 
spectively. Quantitative agreement between experiment and 
theory demonstrates that the resulting repulsive part of the 
depletion forces cannot be interpreted entirely in terms of en- 
tropic arguments but that particularly at small distances (< 
100 nm) attractive dispersion forces have to be taken into 
account. 
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The stability of systems like polydisperse mixtures of 
colloids, natural rubbers, micelles or polymer coils is 
known to be strongly influenced by depletion forces. In 
addition, such forces may also play an important role 
in biological systems, e.g., in promoting the aggrega- 
tion of red blood cells . The understanding of these 
forces being responsible for spontaneous structure forma- 
tion, phase separation or flocculation in such systems is 
highly demanding both from the experimental and the 
theoretical point of view. A first approximate explana- 
tion of depletion forces was given by Asakura and Oo- 
sawa (AO), who recognized that the presence of small 
hard spheres (referred to as macromolecules in the fol- 
lowing) can mediate effective forces between two larger 
objects if their distance is sufficiently small ||. This 
can be easily understood by considering, e.g., a hard 
sphere of radius R suspended in a hard-sphere fluid of 
macromolecules with radius r and bulk number density 
pi, in front of a wall (or another big sphere) at distance 
z, measured between the surfaces of the wall and the 
spheres. For z < 2r, the macromolecules are expelled 
from a region of excluded volume overlap, i.e., their den- 
sity is depleted at that hemisphere of the large parti- 
cle facing the wall compared to the opposite side which 
faces the bulk liquid. Such density gradients give rise 
to an effective osmotic pressure, causing the big sphere 
to be pushed towards the wall thereby allowing the en- 
tropy of the macromolecules to increase. Accordingly, 
these depletion forces result from an asymmetric den- 
sity distribution of the macromolecules around the big 
sphere. Quantitative calculations of depletion forces re- 
quire the knowledge of the density distribution p(r) of the 
macromolecules. Within the crudest approximation p(r) 



is considered to be constant (corresponding to the above 
mentioned AO-approximation), so that depletion forces 
can be easily calculated in terms of excluded volume ar- 
guments and are predicted to be entirely attractive for 
z < 2r. For sufficiently low macromolecule densities this 
approximation is in agreement with experimental results 
If, however, at high concentrations, the structural 
correlation effects in the macromolecular liquid become 
important, p(r) displays an oscillatory behavior at small 
values of z generating repulsive parts of the depletion in- 
teractions as observed recently ||. By means of virial 
coefficient expansions within the Derjaguin approxima- 
tion and by other thechniques in || it has been found 
that the repulsive interaction strength strongly increases 
for large values of R/r. 

So far the interpretation of the corresponding experi- 
ments has been based on assuming hard sphere macro- 
molecule and hard wall interactions which is, however, 
not always justified under experimental conditions. This 
is most evident in the case of highly charged macro- 
molecules where even at volume fractions p{,47rr 3 /3 < 
0.02 electrostatic interactions lead to long-range repul- 
sive and short-range attractive parts in the depletion po- 
tential between a sphere and the wall ^,0. But even 
for systems of uncharged macromolecules or in the high 
salt limit, the model of hard spheres and a hard wall is 
not applicable within the range of attractive dispersion 
forces. 

Here we present measurements of the depletion poten- 
tial of a single colloidal polystyrene (PS) sphere in front 
of a flat glass surface in the presence of non-ionic macro- 
molecules for various values of the macromolecule density 
pi,. In addition we have calculated the depletion potential 
within DFT for a hard-sphere system in which the dis- 
persion and electrostatic interactions between the macro- 
molecules and the wall have been taken into account. We 
find clear evidence for the occurrence of attractive and 
repulsive interaction energies both in our experiment and 
in our calculations. The latter demonstrate that disper- 
sion forces considerably modify the depletion potential 
below 100 nm (~ r). This substantiates the hypothesis 
ifllf that in general for depletion forces both entropic and 
energetic contributions are important. 

We have used TIRM Jl^] to measure the potential en- 
ergy of a colloidal sphere close to a glass surface. In 
this technique, a laser beam is reflected from a solid-fluid 
interface slightly above the critical angle of total reflec- 
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tion thus generating an evanescent wave in the fluid. A 
colloidal sphere located sufficiently close to the interface 
scatters the evanescent wave with an intensity which de- 
cays exponentially with increasing particle-interface dis- 
tance. Thus the scattered intensity, which fluctuates ow- 
ing to Brownian motion, determines sensitively and in- 
stantaneously the separation distance z. In order to ob- 
tain the spatial dependence of the potential energy of the 
particle one has to measure the separation distances sam- 
pled by the colloidal sphere for a statistically long period 
of time. From this, the probability of finding the parti- 
cle at any separation distance can be calculated which is 
related to the potential energy via the Boltzmann distri- 
bution. If the penetration depth of the evanescent wave is 
chosen to be of the order of several hundred nanometers, 
this technique allows one to determine precisely effective 
particle-wall interaction potentials close to an interface. 

Monodisperse PS particles of radius i?=5 pm were sta- 
bilized with sulfate surface groups which dissociate in 
contact with water and thus cause the particles to be 
negatively charged. Owing to the large weight of the PS 
spheres, it was necessary to increase their buoyancy to 
facilitate large enough thermally driven vertical excur- 
sions of the particles during the measurements. There- 
fore the particles were suspended in a mixture of deuter- 
ated water and normal water which caused a reduction 
of the weight by a factor of 8.5. As macromolecules we 
chose poly(ethylene oxide) (PEO) which had been suc- 
cessfully employed in other studies of depletion effects 
before [pp|jl^]. PEO is a non- ionic water-soluble poly- 
mer which forms polymer coils below a critical density 
p c . For the molecular weight which was used in our ex- 
periments (Mw = 10 6 ) the radius of gyration tq has 
been determined by means of static light scattering to 
be r<3=67.7 nm |14| . The maximum polymer number 
density used in our experiments was 23.4 pm~ 3 which is 
smaller than p c by more than a factor of 30. Accordingly 
in the present context the polymer coils can be consid- 
ered as spherical macromolecules to be modeled in our 
calculations as hard spheres {pb <C p c ) with an effective 
radius r=107 nm |fl5j| . 

The sample cell was composed of two silica flats sep- 
arated by a 1 mm thick O-ring and was connected to 
a closed circuit system containing an ion exchange resin 
and a conductivity probe. This setup allowed us to con- 
trol precisely the ionic strength of the suspension. The 
evanescent field was generated by a HeNe-laser beam (A 
= 633 nm, P — 5 mW) which was coupled to the bottom 
plate of the cell under total internal reflection conditions 
by means of a glass prism. The penetration depth was 
chosen to be 400 nm. Highly diluted suspensions with PS 
number densities below 0.5 mm -3 have been used to en- 
sure that only a single PS particle was in the field of view 
and contributed to the scattering signal. The scattered 
intensity of the PS sphere was collected by a microscope 
objective and focused onto a photomultiplier tube. In or- 



der to obtain sufficient data and to minimize statistical 
errors the scattered intensity of a particle was sampled 
over at least 1000 s at a sampling rate of 50 Hz. 

In the absence of depletion forces, i.e., without adding 
macromolecules, the potential of a negatively charged PS 
particle above a glass plate is composed of two parts. To- 
wards smaller distances the potential increases exponen- 
tially owing to the electrostatic interaction between the 
sphere and the fused silica plate, which is also negatively 
charged when in contact with water. For larger distances 
the potential is completely dominated by gravity. It has 
been shown theoretically and experimentally that under 
these conditions the external potential is given by |12 



0*™{z)=Bexp(-Kz)+l3Gz 



(1) 



where B depends on the sphere and the glass plate sur- 
face potentials, k is the inverse Debye screening length 
related to the ionic strength of the suspension, and 
(3 = (fcflT) -1 . G = (pps — Pw)Vg is the weight of the 
particle of volume V suspended in water, pp$ and pw 
are the mass densities of PS and water, respectively, and 
g is the acceleration of gravity. The values for V and pps 
were provided by the manufacturer. 

This potential is shown in Fig. |l](a) as measured (□). 
The PS sphere probes a range of about 300 nm from 
the surface during our measuring time. Due to the large 
weight of the PS sphere with R = 5 /im, we are very sen- 
sitive to changes in the sphere-wall potential below z = 
150 nm; for smaller spheres this potential is dominated by 
the electrostatic repulsion ||] . The solid line in Fig. |l](a) 
is the fit curve according to Eq. (|l|) with the two fit pa- 
rameters B = 4.81 and k" 1 = 17 nm (for T=298 K), the 
latter being in agreement with the electrical conductivity 
measurement of the colloidal circuit. This demonstrates 
that in Eq. (|l|) attractive dispersion forces are negligible. 

However, for the macromolecules dispersion and elec- 
trostatic interactions between them and the glass plate 
have to be taken into account. The latter originate from 
the fact that the charged glass plate induces a dipole mo- 
ment in the neutral macromolecules which interacts with 
the electric field resulting in a repulsion due to the much 
smaller static index of refraction e\ of PEO than that of 
water (e = 81). This leads to 



PEO 



(z) = C exp(-2n(z + r))) - A 



z + r 



(2) 



with the Hamaker constant A = 5 x 10 J |16 and 



C = 12n Ur coshUr) - smh(Kr)) — (3) 

e\ + 2e n A 

with £i w 5 jl7]], the surface electric potential of the glass 
plate Eq 50 mV |l8|], and Sq the permittivity of vac- 
uum. Estimates indicate that for the densities pb con- 
sidered here the attractive dispersion forces among the 
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macromolecules and between them and the PS spheres 
are not important. The agreement, which we find be- 
tween our experimental data and our calculations with- 
out considering the latter forces, supports this approxi- 
mation a posteriori. 

When PEO is added to the system the presence of 
depletion forces generates modified effective potentials 
acting on the PS sphere. Figure [l] shows the mea- 
sured potentials (symbols) <&f Q s t (z) = $%? t (z) + &de P (z) 
when different amounts of PEO are added to the suspen- 
sion. Assuming a homogeneous distribution of PEO over 
the whole circuit, from the experimentally determined 
weights, the PEO number densities are obtained as (b) 
p e b xp = 4.1 /jm" 3 , (c) 6.3 /im" 3 , (d) 8.7 ^m" 3 , (e) 16.4 
/im~ 3 , and (f) 23.4 /im" 3 , respectively. With increasing 
PEO concentration the potential deepens and the posi- 
tion of the minimum considerably shifts to smaller values 
of z in agreement with previous measurements [B . Addi- 
tionally, however, we find a bump in the potential around 
z = 150 nm which increases with increasing PEO densi- 
ties. In order to highlight this feature we have subtracted 
the electrostatic and gravitational contribution (Eq. (Q)) 
from Fig. |l](c) and present this difference in the inset of 
Fig. [I] as measured (A) and as calculated (solid line). 
Within <&dep the bump turns into a pronounced max- 
imum with a height of approximately 0.75 fcgT which 
gives rise to a repulsive depletion force of the PS sphere 
upon approaching the wall before it finally becomes at- 
tractive for z smaller than 150 nm. 

Our theoretical calculations for $>f f (z) are based on 
a novel and versatile DFT approach JnJ. It has been 
shown that in the limit of present interest, i.e., when the 
density of PS goes to zero, the depletion potential $>dep(z) 
can be determined in two steps solely from the density 
profile of the macromolecules Ppeo(z) without PS sphere 
and from the geometrical shape of the PS particle, i.e., 
its radius R [|19|. First, the density profiles of the poly- 
mers ppeo( z ) subjected to the external potential given 
by Eq. (^|) have been obtained by freely minimizing the 
Rosenfeld functional PQ] . With all parameters in the ex- 
ternal potential and the effective hard sphere radius r of 
the macromolecules fixed, the only adjustable parameter 
in the calculation is the PEO bulk density. In Fig. || we 
show a PEO density profile for pfc=4.1 p,m~ 3 (solid line). 
The density profile reflects the short ranged electrostatic 
repulsion at very small distances and the dispersion at- 
traction leading to a pronounced maximum more than 
50% above the bulk value. The importance of for 
bringing about this structure is underscored by the dot- 
ted line in Fig. ^ which corresponds to = 0. Sec- 
ond, from such PEO density profiles we have calculated 
the effective potentials shown as solid lines in Fig. [l](b)- 
(f). As mentioned above, the only adjustable parameter 
in the calculations is pb which has been varied until best 
agreement with the experimental data has been obtained 
leading to the theoretical values (b) p l h h =3.9 /im~ 3 , (c) 



4.1 [Mm' 3 , (d) 4.5 ^m~ 3 , (e) 11.3 /urn" 3 , and (f) 15.5 
/im -3 , respectively. We find very good agreement with 
the experimental potential data and in particular the 
emerging barrier for increasing pb is reproduced. The 
values of pt agree within a factor of about 1.5 with p? xp '. 
A possible reason for these differences could be due to 
lateral density gradients within the experimental circuit 
not captured by the aforementioned global determination 
oip e b xp . 

We emphasize that the potential barriers in Fig. |l| can- 
not be explained in terms of a hard-sphere and hard-wall 
model, i.e., by purely entropic arguments. This conclu- 
sion is demonstrated in Fig. |^ by the discrepancy between 
the actual depletion potential and the one close to a hard 
wall. 

Another important conclusion which can be drawn 
from our calculations is that, although the dispersion at- 
traction seems to be mainly responsible for the potential 
barrier, the soft electrostatic repulsion plays an impor- 
tant role as well. By changing the salt concentration and 
hence the inverse screening length, one can experimen- 
tally tune the range of the electrostatic repulsion and 
thus modify the potential barrier. This explains why in 
an earlier experiment with a lower salt concentration no 
repulsive depletion forces were found js| . 

So far, we have demonstrated that DFT allows one to 
calculate accurately the depletion force on a suspended 
sphere from the density profile of the dissolved macro- 
molecules in the absence of the big sphere. For small 
macromolecule densities this relationship can be inverted: 
measuring the depletion potential of a large sphere may 
be used to determine the undisturbed macromolecule den- 
sity profile. This is a remarkable result, because the big 
sphere considerably modifies the density distribution of 
macromolecules. Thus depletion force measurements can 
be used as a tool for determining undisturbed solvent 
properties. 

In summary, we have studied the depletion potential 
of a single polystyrene sphere in a fluid of non-ionic 
macromolecules in front of a wall. The potential ex- 
hibits both attractive as well as repulsive parts which are 
not of purely entropic origin but are largely dominated 
by the dispersion forces between the uncharged macro- 
molecules and the wall. These results are in good quan- 
titative agreement with density functional theory calcu- 
lations carried out specificly for this system. 

C.B., D.R., and P.L acknowledge financial support 
from the Deutsche Forschungsgemeinschaft through SFB 
513. 
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FIG. 1. Comparison of the effective potential <J?^,f (z) = <J>^(z) + <E>de P (z) acting on the big PS particle between theory 
(solid lines) and experiment (symbols) for various number densities of the macromolecules: (a) p £ b xp = 0, (b) 4.1 /im~ 3 , (c) 6.3 
~ 3 , (d) 8.7 iinr 3 , (e) 16.4 ^m" 3 , and (f) 23.4 fim . For reasons of clarity the curves have been separated from each other 
vertically by 4 units. 
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FIG. 2. The number density profile Ppeo(z) of the macromolecules in the absence of a PS sphere for the wall potential 
&eM°( z ) given by Eq. |i| (solid line) and for a hard wall (dotted line); the theoretical value pl h = 4.1 pmT 4 corresponds to 
Fig. 0(c) and P peo(z = 0)(2r) 3 = 0.043. 
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FIG. 3. The depletion potentials <&dep(z) between the wall and a big PS sphere calculated from the corresponding PEO 
density profiles ppeo(z) shown in Fig|| A denote the corresponding experimental data. 
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